Enhancement of pairwise entanglement via symmetry breaking 
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We study the effect of symmetry breaking in a quantum phase transition on pairwise entanglement 
in spin- 1/2 models. We give a set of conditions on correlation functions a model has to meet in order 
to keep the pairwise entanglement unchanged by a parity symmetry breaking. It turns out that all 
mean- field solvable models do meet this requirement, whereas the presence of strong correlations 
leads to a violation of this condition. This results in an order-induced enhancement of entanglement, 
and we report on two examples where this takes place. 



Entanglement is a non-locality inherent to quantum 
mechanics and an important resource for quantum op- 
tics and quantum information processing. It also has 
attained a lot of interest in the last decade from con- 
densed matter physicists, especially since a connection 
between quantum non-locality and quantum phase tran- 
sitions was proposed P, [3] . This initiated a vast analysis 
of quantum critical models with respect to their entangle- 
ment features based on the few computable entanglement 
measures at hand 0, i, i, S, 0, i, i, E El, Elll, \M, 



15l . Il6l . |17[. One of the entanglement measures probed 
on such models is the concurrence [l^, [3, [lOl , a measure 
of pairwise entanglement for spins- 1/2 (qubits). Most of 
the results obtained in the literature make use of certain 
symmetries of the model Hamiltonian which reflect in the 
reduced density matrix. In presence of quantum phase 
transitions, a symmetry is typically broken in the ordered 
region with non- vanishing order parameter. The effect of 
the broken symmetry on the entanglement largely has not 
been taken into account in the existing literature (but see 
e.g. for the local entropy). It was exactly this subject 
which has been addressed by Syljuasen in [2l|, [13] for the 
quantum-Ising and the ferromagnetic XXZ model, where 
a parity symmetry is broken at the transition. He pro- 
vided a condition on correlation functions that guaran- 
tees the invariance of the concurrence under symmetry 
breaking. The models under consideration there fulfill 
this condition. The focus of the present work is to inves- 
tigate the regime where Syljuasen's condition is violated 
as it happens for the XY model in transverse field. 

In this letter, we first summarize the main result of 
Rcf. [2l|. Then, we derive the condition for two-point 
correlation functions and form factors that guarantees 
the invariance of the concurrence under parity symmetry 
breaking where the Syljuasen condition docs not apply. 
Finally, we illustrate our result, showing that mean- field 
solvable models satisfy this condition and demonstrating 
that symmetry breaking in general must be given account 
for, both for the XY spin chain in transverse magnetic 
field and the Lipkin-Meshkov-Ghck (LMG) model. 



We will next review the main result of Ref. l21|, adopt- 
ing for the sake of cross-readability the notation used 
in this work. The model Hamiltonians H under con- 
sideration have a parity symmetry or global phase flip 
symmetry, meaning that the eigenstates can be chosen 
as superpositions from states all having the same par- 
ity of flipped spins (i.e.: odd or even). It is reflected by 
[Hi erf, i/] = 0. The general 2 sites density matrix of 
such a system is given by 



( A a a F\ 

a B C b 

a C B b 

\F b b D J 



(1) 



where indices for the entries of the density matrix have 
been omitted. The entries of p are related to spin corre- 
lators as follows: a = {(S'=) + 2 {S'^S'))/2, b = ((5^) - 
2(S'^S'"))/2, ^= {S'S') + {S') + 1/A,B = l/4-{S'S'), 
C = (S'^S'^) + {Sysy), D = {S'S') - {S') + 1/4 and 
F = (S'^S'^) - {SySy). The symmetry-breaking man- 
ifests itself in a, 6 7^ 0. For a = b = Q. the square 
root of the eigenvalues of pp are B ± C and VAD ± F. 
The concurrence of a 2-site reduced density matrix p is 
computed from the positive semidefinitc matrix R := 
pi(7y (g) ay)p*{ay (g) tTy ) as C = max{0, 2X,nax - tr VR}, 
where Xmax is the maximum eigenvalue of \/R. Thus the 
concurrence is C = 2max|o, |C| - VAD, \F\ - =: 

2max{0,Ca/,C/}. We refer to the concurrence as be- 
longing to the anti-ferro- or ferromagnetic sector when- 
ever the largest eigenvalue of ^/R is from this sector and 
denote it by Caf and Cf, respectively. Throughout all the 



letter will we adopt the notation {S"'S'^) := (^S°'S'^j and 

{S") := {{Sf) + {S^))/2. For the models considered here, 
the Hamiltonian is real and so is the density matrix. 

In order to study the effect of symmetry breaking, 
accompanied by non-zero {S^S^} and order parameter 
(S^), one has to extract the dependence of the eigen- 
values of i? on a and b. This task simplifies since the 
eigenvector (0,1,-1,0) (the singlet Bell state) and its 
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eigenvalue B — C are unchanged. If C > 0, this robust 
eigenvalue can not become the largest one. 

The exact factoring of a cubic polynomial is not very 
handy and has been cleverly avoided by Syljuasen, when 
concentrating on the coefficients of the characteristic 
polynomial — g2X^ + giX — go of R with go = (xyz)^, 
gi = cc^y^ + .T^z^ + y^z^, and g2 = x"^ + + z"^ rather 
than studying the eigenvalues of y/R themselves. The 
connection between the gis and the entries of the re- 
duced two-spin density matrix is 21 1 



go ^ {a'^ - 4'yS)P - Afii^a - i^i'^S 
gi = a'^ + 2aP - Ai^iv - A-^d 
52 = 2a -1-/3 



(2) 
(3) 
(4) 



where a = F"^ + AD - 2ah, (3 = [B + Cf - 4a5, 7 = 
DF-b^, S = AF -a^, 11 = aD-b{B + C - F) and v = 
a{B+C—F)—bA. The concurrence would then simply be 
expressed as C = max{0, k— (i? — C)}, with k z — x — y 
and z the largest eigenvalue of ^/R. The condition that 
the concurrence be unaffected by symmetry breaking is 
then [Hi 



2kV5o = - gi 



(5) 



Interestingly, for dominating Cj, i.e. when vAD -|- |F| is 
largest, this condition is identically satisfied which is the 
main outcome of Ref. 
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We will now extend the analysis of Rcf. [2H to the 
regime where Caf dominates, but still C > 0. The case 
C < and B — C being the largest eigenvalue of ^fR^ is 
treated the same way, just replacing nhy k := x + y + z. 

It is a generic feature of quantum spin models (with 
competing interactions) that there is a point in cou- 
pling parameter space at which the ground state factor- 



15l l23l . I24| . At such a factorizing field any mea- 



izes 

sure of entanglement must vanish. If at both sides of a 
factorizing point the concurrence is non-zero, there are 
two scenarii: either the concurrence is smooth at hf or 
the sector giving support to the concurrence necessar- 
ily has to change at this point. It is the latter sce- 
nario which actually takes place in the models exam- 
ined so far in the literature as e.g . the transverse XY 
model [l|,[2i], the XYZ model [110, and the Lipkin- 
Meshkov-Ghck model (LMG) [sl. l26t"~the Syljuasen con- 
dition applies between the factorizing and the quan- 
tum critical field. The explicit form of the concurrence 
for parity symmetric models and translational symme- 
try (which includes all above cited studies) implies that 
Caf = C/ = at the factorizing field, and they cross \xq.- 
less2(5'2^S'2')-((5'^5'^)-(S'^)2)/(2(5^)2) - [h-hfY with 
even n > (for symmetry breaking field in x-direction) . 
If only one crossing occurs, then the concurrence is robust 
against symmetry breaking between critical and fac- 
torizing point. Below the factorizing field Caf dominates 



and the concurrence will in general be affected by the 
broken symmetry. Nonetheless, to the best of our knowl- 
edge no such case has been reported in the literature, 
which motivates the search for further conditions that 
make the concurrence robust against symmetry breaking 
or for examples where the concurrence is affected. We 
will show in this letter that both cases occur. 

Having taken the square of Eq.® and inserted n 
B + C — 2\J AD, the resulting expression is conveniently 
expressed in the new variables A := b/a and a, leading to 

(6) 



32a^{X ~ Xof [1 + a^(A - Ai)(A - A2)] = 



where Ao = \/D/A, Ai + A2 = 2k/A and (Ai - A2)^ = 
4k{B + C)/A\ 

Both Ai and A2 are real for k > 0, which is mandatory 
for a non-vanishing concurrence. The real solutions to 
Eq. ([6]) arc constraints on the correlation functions of 
the model which insure the concurrence to be unaffected 
by symmetry breaking. The solution A = Ao can then be 
recast into a simple condition: 



{1 + 4 (S^S^))^- 16 {S^f _ - 2 (S-S^) 



1 



-4(5^) +4(5'^ S*^) (5^) + 2 (S'-^S'^) ■ ^'^^ 

Note that this condition is automatically satisfied for 
mean-field solvable models as a direct consequence of the 
factorization of the two-point functions in this case. 
Another real solution to Eq.® exists if a < —1 and 



Ai — A2 
X1+X2 + 2 



1 

±- 

TT 



(8) 



to insure the reality of both a and b. This leads to the 
following transcendental equation 



B + C , 

"(;^W^^i"' 



1 



A{k + A) 1 
A{k + A) + 2k? 7r|a| 



(9) 



Eq. Q has a unique real solution for arbitrary A, B 
and K. The rewriting of this condition in terms of corre- 
lation function is lengthy but straightforward. There is 
numerical evidence that a < — 1 only occurs for a two-site 
reduced density matrix of rank 3 or less. Interestingly, 
rank 3 and 2 are assumed for the (2-magnon) W state 
T.iW^JStS+ and GHZ (cat) state \^) + |tr) of the 
whole chain, where Wij € C and Iff) and |-|J-) are in op- 
posite direction fully polarized states. If Eqs. ([7]) and 
([H])-® are violated, the concurrence will be affected by 
symmetry-breaking; as condition (|5l) is then violated. 

LMG model The LMG model [27| in transverse field 
is described by the Hamiltonian 



where the cr^ are the Pauli matrices for the i-th lattice 
site and N is the chain length. Any two spins are in- 
teracting with the same coupling strength. The pref- 
actor 1/N leads to a finite free energy per spin in the 
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thermodynamic limit. H commutes with the total spin 
and the spin-flip operator Yii o'l for any anisotropy pa- 
rameter 7. This discrete symmetry is broken at the 
quantum critical point, he = 1. The factorizing field of 
this model is /i/ = 1^/7. Due to its infinite-range cou- 
pling, the mean-field approximation for the LMG model 
becomes exact in the thermodynamic limit psj . In this 
limit, condition ([7| is satisfied for arbitrary transverse 
field. The non-trivial behavior of the entanglement of 
the LMG lies in the finite-size corrections to the mean- 
field solution SEi]. 



Hamiltonian |22 



symmetric (even) 

symmetry broken (even + odd) 

AC„ (xlO) 




FIG. 1: Rescaled concurrence for the LMG model: The even 
sector (full line) and the symmetry broken (dashed line) con- 
currence is shown for TV = 100 and 7 = 0.5; the inset shows 
the difference between the rescaled concurrence in the Z2- 
symmetric case and in the Z2-broken case for fixed magnetic 
field h — 0.2 < hf. It scales as ao/N + ai, with ao = 0.651622 
and tti ~ 10~^. The dash-dotted curve is this difference 
rescaled by a factor of 10 as a function of the magnetic field. It 
goes to zero at /i ~ 0.5. At such a point, one of the conditions 
for the correlation functions of the model must be satisfied. 



In Fig. [T] we report numerical data for the rescaled 
concurrence, Cr = {N — 1)C, for the even (black full 
curve) and K2-broken (blue dashed curve) ground state. 
In the thermodynamic limit {N goes to 00 while the free 
energy per spin remains finite), both curves coincide. The 
analytic expression of this limiting curve in the region 
< h < hf has been obtained in [26|: 



lim C 

AT— foo 



Sym 



1 



■7 



1 - h^ 



(11) 



For finite N, we clearly see a deviation of the concurrence 
calculated with a state which violates the spin-flip 
invariance from the concurrence evaluated with a state 
which conserves the latter symmetry. This is indeed a 
finite-size effect which is beyond the scope of the first 
order quantum corrections performed in Ref. [2^. In 
fact, the difference between C'^*'™ and 0^^°'^'''^ is of order 
1 /N"^ as shown in the inset of Fig. [T] 31 1 . 

XY model The transverse XY model is defined by the 



HxY — — - 



where h is the transverse field strength and 7 the 
anisotropy parameter. Below the quantum critical field 
/ic = (1 + 7)/4 the parity symmetry is broken; the fac- 
torizing field is at hf = Wc numerically com- 
puted the ground state nearest neighbor (n.n.) concur- 
rence by means of the Density Matrix Rcnormalization 
Group [29|. This powerful numerical technique finds an 
optimal truncated bases of size m ^ 100 — 200 ^ 2^ 
to describe the spin chain wave function keeping the de- 
sired precision |30[ • In Fig. [2] wc show the results of 
the numerical simulations for the concurrence of the cen- 
tral spins (in order to minimize boundary effects) of a 
chain of length N = 199 and different anisotropy values. 
The concurrence for the even sector C^v™- and broken 
symmetry ground state C-^rofcen then compared for 
different anisotropy values. The concurrence vanishes at 
the factorizing field and its derivative with respect to 
the transverse field at the critical point diverges as ex- 
pected [l[. The even and odd sector concurrences are 
found to coincide for the parameter range we considered 
(data not shown). We used an additional field ft-s along 
Ux to break the symmetry; our results are stable with re- 
spect to changes in the field strength hs G [10"'* : 10~*] 
and the field direction in the XY plane (data not shown) . 
Below the factorizing field the concurrence of the bro- 
ken symmetry and even ground state are clearly differ- 
ent while they are equal between the factorizing and the 
critical field as there the condition ([5]) holds [2^. In 




FIG. 2: N.n. concurrence for the anisotropic XY model for 
N = 199, 7 — 0.3, 0.5, 0.7 (circles, squares and diamonds 
respectively), m = 200 in the even sector (full) and for the 
state with hs ~ 10~® (empty). 



Fig. [3] we plot the difference between the broken symme- 
try concurrence and the concurrence in the even sector 
from the data of Fig. [51 As expected, for h/hf > 1 the 
difference is zero while on the left of the factorizing field 
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{h/hf < 1) the difference is not negligible. The inset of 
Fig. E] reports the finite size scaling of {C'^y"^ - C^''°'='="i. 
Differently from the LMG Model, the broken symmetry 
effect on the concurrence does not vanish in the ther- 
modynamic limit. It approaches a non-zero value with 
oscillating behavior. This means that for the XY model 
conditions ((Zl)-® and thus Eq. ([5]) are violated below hf. 
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FIG. 3: Difference of the n.n. concurrence in the even sector 
and for the symmetry broken ground state from the data of 
Fig. [2] as a function of the transverse field. The maximum 
relative deviation amounts to around 10%, decreasing with 7 
for sufficiently small h. Inset: Finite size scaling for 7 — 0.7 
and h/hf — 0.8 (diamonds) and limiting value (dashed line). 



Conclusions We have found conditions on spin cor- 
relation functions, which ensure that the concurrence is 
invariant respect to breaking of a parity symmetry. They 
are necessary and sufficient in a regime complementary 
to where the relations from Ref. 2l| do apply. For mean 
field exact models, one of the conditions (Eq. [7]) is sat- 
isfied and hence will the rescaled concurrence be unaf- 
fected by parity symmetry breaking. The further condi- 
tions (Eqs. in [9]) only emerged for the reduced two-site 
density matrix having rank three or less. Interestingly, 
this occurs for the whole system being in a W-type state 
with two running flipped spins in the ferromagnetically 
polarized state and a W or GHZ-state, respectively. It 
would be interesting to verify, whether such states may 
satisfy the corresponding condition on correlation func- 
tions; W-states would lead to a long range concurrence, 
which has been observed in [l^ close to the factorizing 
point. Numerical studies reveal that the concurrence is 
affected by symmetry breaking for the LMG and the XY 
model. As a consequence, the conditions (III)-® on the 
spin correlations are violated up to perhaps a single value 
for the magnetic field in the LMG model. We gave cer- 
tain conditions that might admit scenarii different from 
a robust concurrence between he and hf. Their investi- 
gation could provide important insight in the interplay 
of entanglement and quantum phase transitions. 
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